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In order to accept very forward angle scattering particles, Jefferson Lab HKS exper- 
iment uses an on-target zero degree dipole magnet. The usual spectrometer optics cal- 
ibration procedure has to be modified due to this on-target field. This paper describes 
a new method to calibrate HKS spectrometer system. The simulation of the calibration 
procedure shows the required resolution can be achieved from initially inaccurate optical 
description. 

1. Introduction 

Jefferson Lab HKS experiment aims at obtaining high resolution hypernuclear spec- 
troscopy (e.g., ~ 400 KeV FWHM for fB) by (e,e'K) reaction (ref.JIj). It utilies CE- 
BAF primary electron beam. The electroproduction provides a advantage over hadronic 
production mechanism because the momentum and position spread of the primary elec- 
tron beam is significantly less than secondary hadron beams. The kaons are detected by 
High Resolution Kaon Spectrometer (HKS) with a momentum resolution of 2 x 10~ 4 . In 
order to maximize hypernuclear yield, the target is placed just before a zero degree dipole 
magnet (Splitter). The field of the splitter is used to seperate scattered e' from K + and 
allow the spectrometer to accept both e' and K + from very forward angles. 

To calculate hypernuclear excitation energies, we need to determine scattered kaon and 
electron momenta and angles. The kaon momentum is the dominant factor for the final 
missing mass resolution. In order to obtain the proposed missing mass resolution, it is 
important to optimize the spectrometer optics to get the required momentum and angle 
reconstruction resolution. This calibration procedure is complicated for HKS spectrometer 
Because of the existence of splitter magnet right after the target. 

Table 1 

required momentum and angle reconstruction resolution 




The spectrometer optics can be conviniently expressed in a matrix formalism (ref.[2]), 
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Usual double arm spectrometer HKS spectrometer system 




Figure 1. Schematic view of magnetic optical properties of normal double arm spec- 
trometer system (left) and HKS spectrometer (right) with an on-target dipole magnet 
(Splitter). 



which is a set of polynimial coefficients related target angles and particle momentum 
to focal plane measurements. Conventionally, the method for calibrate matrix for angle 
reconstruction is by using a special designed collimator called sieve slit with a array 
of holes drilled onto metal template (ref.[2]). The sieve slit is mounted in front of the 
spectrometer entrance so all the particles passing through the holes will have defined 
incident angles. 

But in the setup of HKS experiment, because target is inside the field of Splitter, the 
sieve slit plate can only be mounted on the Enge and HKS dipole entrance after the 
Splitter (figHJ). The particle trajectories are already bended by the splitter field before 
they pass through the sieve slit holes. Thus a hole of sieve slit no longer select a single, 
fixed angle. The usual way of angle reconstruction calibration has to be modified. 

For momentum calibration, because Enge, HKS two arms and the beam dump line 
are all coupled through splitter field, thus suitable kinematics setting for typical single 
arm momentum calibration (e.g., by production of carbon excited states) is difficult to 
implement. 

In this paper a high precision optical calibration method for HKS spectrometer is de- 
scribed. The procedure can be applied to optical calibration of other magnet spectrometer 
system with on-target magnetic field. The optical calibration procedure will use the data 
from Enge and HKS sieve slit runs for angle reconstruction calibration and CH 2 target 
runs for momentum calibration. Before the experiment, the procedure has been tested by 
simulation with initial optics. The simulation is carried out by using measured field map 
of the spectrometer magnets and tracking particles through the field by program GEANT 
(ref. |Hj). The calibration procedure for angle, momentum, kinematics reconstruction 
and beam raster correction will be described from Sec.l to Sec. 5 seperately. The iteration 
of these procedure is described in Sec. 6. The simulation and result of the calibration 
procedure from is presented in sec. 7, followed by a summary in section 8. 
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2. Sieve Slit Calibration 

In HKS experiment, although the Splitter field smears the one to one correspondence 
between target angles and sieve slit hole position. The trajectory of a particle which 
passing through a given sieve slit hole and scattered from a point target is a functional of 
the paticle momentum and splitter optics. Thus with known splitter optics, the scattering 
angle can be expressed as a function of particle momentum and sieve slit hole position. 
The optics of Splitter can be obtained by measuring the field distribution of the magnet 
and tracking particles through the field by programs such as GEANT (ref.jSJ). From the 
point of view of total integral of Bdl, the Splitter contribution is only ~ 8%, because the 
Splitter is just used to seperate particles. The contribution of field measurement error is 
negligible for target angle reconstruction. 

Based on this observation, the sieve slit(SS) calibration procedure is: 

1. Fit function F s2 t- target angles as a function of SS hole positions and particle 
momentum based on Splitter optics. 

2. Separate events of the SS calibration run data hole by hole, thus determine corre- 
sponding SS hole center position for each event. 

3. Using function F S 2t, calculate corrected target angles (xptar, yptar) for each event 
from SS center position and momentum. 

4. Fit transfer matrix from focal plane to target: {xfp, xpfp, yfp, ypfp) — >(xptar, 
yptar). 

A simulation of the SS calibration procedure with initial inaccurate Splitter, Enge optics 
reaches final resolution for e' (Enge) target angles (a): xptar: 1.8 mr, yptar: 1.1 mr, for 
kaon (HKS) angles: xptar: 0.4 mr, yptar: 2.4 mr. 

3. Momentum Calibration 

The method of HKS momentum calibration is making use of the known masses of A,S° 
hyperons and the narrow width of \ 2 B hypernuclear ground state (ref.jH). The hyperons 
is produced from hydrogens in the CH2 foil (5 mm in thickness) The momentum recon- 
struction matrices for Enge and HKS arms are optimized simultaneously by minimizing 
the Chisquare defined as the sum of squared mass differences between the calculated mass 
and the PDB values. A Nonlinear Least Square method is used to optimize the matrices 
for electron arm and kaon arm simutaneously. 

Although there is emulsion measurement of ^ 2 B GS binding energy at 11.37 ± 0.06 
MeV(ref. jS]), the associate statistical error and the composite nature (2~/l~) of the state 
may cause the peak center to deviate from this value (ref. [H]). So in the calibration 
process, the peak position is a adjustable parameter and the calibraion is carried out in a 
iterative way, the GS peak center is adjusted after each iteration according to the actual 
fitted position of calculated \ 2 B GS distribution, until an minimum Chisquare is reached. 

The procedure is: 
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1. Calculate missing mass using the existing angle matrix and the initial (un-calibrated) 
momentum matrix. From CH 2 target runs, select A,E° events, from Carbon target 
runs, select \ 2 B GS events from a window located around the center of the missing 
mass peaks. 

Inside the 5mm thick CH2 target, the energy loss and radiative processes are sig- 
nificant. It shifts the missing mass distributions and forms a tail. The width of 
the event window has to be selected appropriately, because a too wide window will 
dilute the sample with too many radiated events while a too narrow window will 
cause a insufficient statistics in the fitting. 

2. Define Chisquare as the sum of squared mass differences Am- between the calculated 
mass and the PDB values or the initial assumed \ 2 B GS binding energy: 

x 2 = J2 W A™% 

where Wi is the relative weight of A,E and \ 2 B GS events. 

3. x 2 is a function of missing mass emiss, thus a composite function of e' and kaon 
momenta 5p e and 5pk- 

x 2 = f(emiss(5p e ,5pk) 

Minimizing x 2 by by a Nonlinear Least Square method to optimize the Enge and 
HKS momentum reconstruction matrices(ref.[7j). 

4. fit the actual peak center of the calculated ^ 2 B GS missing mass distribution. Using 
this value as the new ^ 2 B GS binding energy. Go back to first step. Iteration until 
an minimum Chisquare is reached. 

4. Kinematics calibration 

The purpose of kinematics calibration is to find the offsets of beam energy, Enge and 
HKS central momentum from nominal values. The method is also to use the known masses 
of A,S° and \ 2 B ground state, same as momentum calibration. The x 2 is defined the same 
way as in momentum calibration. But for kinematics calibration, the % 2 is a function of 
beam energy offset AE e0 , e' central momentum offset Ap e > and kaon momentum offset 
Apk, corresponding to the zero order terms in the momentum matrices. At correct offsets, 
this x 2 w iH De minimized. 

If we define sum of kinematics offsets as: 

Ap kin = AE e0 - Ap e , - Ap k 

The x 2 w iH dominantly depend on Aphin- Correct offsets can be found by scanning on 
Apkin over all possible values to locate the minimum Chisquare. 
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Figure 2. Comparison of Enge angle (left) and momentum(right) reconstruction resolution 
with rastered beam,with and without raster correction. Both resolutions are shown when 
the correction function g fitted as a 3-order and 4-order polynomials. The rightmost point 
is the resolution after raster correction with wrong Enge optics. 

5. Raster Correction 

To avoid burning the CH2 target, beam will be rastered 5mm x 5mm on the target 
by a pair of bending magnets. If the raster effect is ignored, a direct reconstruction 
with matrices obtained for point target gives much worse resolution, for example, the 
5p resolution for HKS is 1.16 x 10 (rms) without raster, it becomes 6.64 x 10" 4 with 
5mm x 5mm raster. 

To do raster correction, the beam positions on target need to be determined by raster 
magnets X and Y current event by event. For rastered beam, The reconstructed target 
quantities can be expressed as functions of focal plane quantities and beam positions at 
target. It can be written as, for example: 

£>P = f(xfp, xpfp, yfp, ypfp) + g(xtar, ytar) 

We seperated the dependence on beam positions from focal plane quantities so that the 
raster correction function g depends only on beam position. It can be fitted from initial 
spectrometer optics. Even if the actual optics of the spectrometer system may deviate 
from the initial optics, its effect on raster correction function g is negligible so we can still 
use the initial function form (figOand figOJ). 

6. Iterative Calibration 

It is obvious that the individual calibration procedure mentioned above actually are 
correlated. For example, the calibration of target angles need the input of particle mo- 
mentum while to extract missing mass, we need to know the scattering angle. Thus the 
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Figure 3. Comparison of HKS momentum reconstruction resolution with rastered 
beam,with and without raster correction. 

individual calibraion procedures should be carried out in an iterative way until conver- 
gence is reached. The general procedure will be: 

1. Separate Enge and HKS sieve slit calibration to obtain angle reconstruction matrices 
using original momentum reconstruction matrix. 

2. Raster correction for CH 2 target data using existing raster correction matrix. 

3. Kinematics calibration to find deviation of beam energy, Enge and HKS central 
momenta from their nominal values. The central angles of the spectrometers are 
not used in the missing mass calculation and they are not needed for fixed-angle 
spectrometer. 

4. Two arm coupled momentum calibration to obtain momentum reconstruction ma- 
trix. 

5. After getting the calibrated momentum matrix, go back to first step again to do 
SS calibration using the new momentum matrix. Iterate until best resolution is 
reached. 

7. Simulation of the Calibration Procedure 

Simulated events are used to test the calibration method. The simulation take account 
of target physics processes, spectrometer optics and detector resolution. 

The simulation of target processes is adapted from SIMC (ref.jBj), a physics simu- 
lation program for Jefferson Lab Hall C basic equipment. The ionization energy loss, 
Bremsstrahlung, multiple scattering and beam spread effects are included in the simula- 
tion. The simulation uses 5mm thick CH 2 target and 0.1 g/cm 2 Carbon target. 



No Raster Correction 



3-order 



4-order 



No Rast 5r 
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Simulated events from target are send through The spectrometers. RAYTRACE (ref.^J) 
is used to simulate events in the electron arm. A measured magnet field map is used to 
track particles through kaon arm by GEANT. Only optical properties (no physical pro- 
cesses) are considered in this step. 

At focal plane of the spectrometers, the positions and angles are smeared according to 
detector resolution and multiple scattering. For Enge, the resolutions (a) are: 86 fj,m (x), 
0.7 mr(xp), 210 /im(y), 2.8 mr(yp). For HKS, they are: 160 fxm (x and y), 0.33 mr (xp 
and yp). 

The simulated focal plane events are then reconstructed back to the target using back- 
ward reconstruction matrix. The beam energy and particle momentum are corrected by 
the average energy loss in the target. The corrections are 0.4923 MeV (beam energy), 

0. 6015 MeV (e' momentum), 0.4890 MeV (kaon momentum) for CH 2 target, 0.0820 MeV 
(beam energy), 0.1026 MeV (e' momentum), 0.0866 MeV (kaon Momentum) for Carbon 
target. 

The simulated A, S° and )?B missing mass spectra are shown in figEJ The predicted 1 
state of \ 2 B at 2.73 MeV is also included in the simulation for comparison. The quasi-free 
carbon background on CH 2 target is simulated with A and E° events by a assumed S/B 
ratio of 6:1 under A peak. The \ 2 B GS has a missing mass resolution of 397 KeV (FWHM) 
with correct optics. 

The simulated Enge focal plane sieve slit correlation patterns are shown in figH3 FigEI 
shows the yfp vs. xfp correlations for each X-column of sieve slit holes. 
The simulated HKS focal plane sieve slit patterns are shown in fig|7| 
To test the calibration method by simulation, We intentionally use wrong spectrometer 
optics in reconstruction. The splitter field is changed from nominal value of 1.546 Tesla to 
1.550 Tesla in forward simulation for e' events, but still using the nominal reconstruction 
matrix in reconstruction. In addition, the beam energy, electron and kaon arm central 
momenta have offsets from nominal values, they are 1.2 MeV, -0.7 MeV and -0.9 MeV 
respectively. The missing mass distributions before the calibration is shown in fig. |HJ 
The missing mass distributions after the calibration is in fig. El The width of the \ 2 B 
GS improved from 1.56 MeV (FWHM) to 0.418 MeV (FWHM) after the calibration, and 
the center of GS peak is within 20 keV of the expected GS mass obtained by emulsion 
experiment. 

8. Summary 

With the use of a dipole magnet on target, the normal spectrometer calibration pro- 
cedure has to be modified. We have developed a new method for calibration of the 
reconstructions for target angle, momentum and kinematics. Starting from a initially in- 
accurate optics, by iteration of these calibration procedures, the spectrometer optics can 
be calibrated to the required resolution. 
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Figure 4. Simulated A, S° and ^ 2 B missing mass spectra with Carbon QF background 
from CH 2 target. 
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Figure 5. Simulated Enge focal plane sieve slit correlation and hole patterns. 
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Figure 6. Simulated Enge focal plane yfp vs. xfp correlations for each X-column of sieve 
slit holes. 
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Figure 7. Simulated HKS focal plane sieve slit correlation and hole patterns. 
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Figure 8. Reconstructed A, S° and \ 2 B missing mass spectra with wrong optics before 
the calibration. 
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Figure 9. Reconstructed A, S° and \ 2 B missing mass spectra with optical matrices cali- 
brated. 
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